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Abstract. We present microscopic calculations of low energy scattering ob- 
servables in all possible four nucleon systems: n- 3 H, p- 3 He and p- 3 H. Re- 
sults were obtained by solving Faddeev-Yakubovski equations in configuration 
space, appropriately modified to include Coulomb and three-nucleon forces. 

1 Introduction 

Three- and four-nucleon systems are the testing ground for studying the nuclear 
interaction. If modern NN potentials have reached a very high degree of accuracy 
in describing the two-nucleon data, they fail already to account for the binding 
energies of the lightest nuclei. The use of three-nucleon forces (3NF) is - with 
some exceptions - mandatory. By adjusting its parameters one can obtain 

a satisfactory description of nuclear binding energies up to A=10 5 . 

Apart the effect of rescaling nuclear thresholds, low energy three-nucleon 
scattering observables are quite insensitive to 3NF. Furthermore, since 3N spec- 
tra does not contain any narrow resonances, three-nucleon dynamics is relatively 
rigid once deuteron and triton binding energies are fixed. The four nucleon con- 
tinuum is a challenging few-body nuclear problem. Its interest lies not only in 
the natural progression that it represents, towards a systematic description of 
nuclear systems with increasing complexity, but also in the richness of the A = 4 
nuclear chart itself (see Fig. Four nucleon problem, we believe, implies a 
qualitative jump in respect to the A = 3 case. 

2 Theoretical tools 

We use the Faddeev-Yakubovski (FY) equations in configuration space to tackle 
the four-particle problem. However in their original form these equations are 
applicable only for systems interacting via short-range pairwise forces. An elegant 
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Figure 1. Four nucleon energy spectra. Single lines indicate particle thresholds, whereas dashed 
lines represent particle decay unstable excited states (resonances). 



way to include 3NF into FY equations was suggested in for four identical 
particles. These equations reads: 



(E-H -V 12 )Hfl = V 12 P[(l + Q)Kt 2<3 + H^ 



(E-H - V l2 ) Kf 2<3 = V l2 {P + + P-) [(1 + Q)Kf 2i3 + Hfl] + V 12 ,3* 

with P + , P~ , P and Q being permutation operators: 

P+ = (P-)- = P 23 P 12 ; Q = eP 3A ; P = P 13 P 24 = P 2A P 13 . (2) 

and 

V = [1 + (1 + P + + P-)] Q(l + P+ + P-)Kf 2 3 + (1 + P+ + P")(l + P)H?t 

the total wave function. 

Equation (^Q) becomes however non appropriate once long range interaction, 
in particular Coulomb, is present. In fact, FY components remain coupled even in 
far asymptotics, thus making the implementation of correct boundary conditions 
hardly possible. In order to circumvent this problem we have split Coulomb 
potential {Vc) into long {V c ) and short (V c sh ^) range parts by means of some 
arbitrary cut-off function similarly as it was done by Merkuriev for 

3-body equations [I]: 

v («)-^*>(«,,,,) + v^(«,,,,)i S, { *' V i~f*'V )V °$l , y 

V c (x, y, z) = [1 - x{x, y, z)\ V c {x) 

(3) 
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x(x,y,z) is a smooth function which equals 1 for small x values or when x <C 
(y, z) and vanishes for x 3> (y, z). Then the 4-body FY equations can be rewritten 
in the form: 

(E-H -V^ h) -J2 V c ] ) K k3 = vlt\P + +P-)[(l + Q)K 4 12 ,3 + H 3 1 l]+V 12 ^ 
iE-H -vW--£v<P)H% = V^P[{l + Q)Ki % s + H%] (4) 

We use for short the following notations: 

V^ h) = (v {N) + V^) i2 . (5) 

In this form FY equations become assymptoticaly uncoupled and appropriate 
boundary conditions can be easily implemented. 

Equations Q are solved by making partial wave decomposition of amplitudes 
Kf 2>3 and Hit. 
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The partial components K.f and Tif ' are expanded in the basis of three- 
dimensional splines. Thus integro-differential equation (j3J) is converted into sys- 
tem of linear equations. For detailed discussion on the equations and method 
used one can refer to HOI. 



3 Results 

3.1 n-^H elastic scattering 

The elastic scattering of neutrons on tritium is the simplest 4N reaction. Treat- 
ment of this system doesn't require calculation of the costly integrals due to 
the Coulomb terms in eq. (@J). Nevertheless this system contains two, spin de- 
generated, narrow resonances at low energy (E cm ~ 3 MeV). Ability of nuclear 
interaction models to describe scattering cross sections in this resonance region 
is still an open question [HJ |HJ - 
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Table 1. 4N scattering lengths calculated using different interaction models. 



Semi-realistic MT I-III potential was shown to be very successful in describ- 
ing total as well as differential cross sections [7] at the energies bellow n-n-d 
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threshold. Situation is less obvious for realistic potentials, which require much 
larger partial wave basis (PWB) to obtain converged results. Recently we have 
managed to considerably enlarge our PWB compared to 8. . Despite of having 
some effect on the negative parity phase shifts, the total cross section near the 
resonance peak (E cm = 3 MeV) has not been improved (see Fig. EJ. We should 
notice, however, that the J n = 2~ phase shifts, the most relevant due to its 
statistical factor, are not fully converged yet. They still show some trend to in- 
crease, however we estimate that it can not suffice to reproduce the experimental 
data. 




Figure 2. Differential cross sections versus centre-of-mass angle for n+ 3 H at E cm = 2.625 MeV 
(figure on the left) and for p+ 3 He at E cm = 4.1325 MeV (figure on the right) 

By including UIX-3NF we have managed to reproduce the experimental zero- 
energy cross sections (see Table 0, which are overestimated by realistic NN 
interaction without 3NF. However their effect near the peak remains very small. 
Other smaller discrepancy exists near the minima of total cross sections (E cm ~ 
0AM eV). In this region, determined by positive parity states, both MT I-III 
and Av.l8+UIX overestimate experimental data. This suggest that low energy 
parameters (0 + and 1 + scattering lengths and effective ranges) are not described 
sufficiently well. 

3.2 p- 3 He elastic scattering 

This system is an isospin partner of n- 3 H. Yet calculations were done uniquely 
with MT I-III model, except of J n = + scattering lengths for which our Av.18 
and Av.l8+UIX predictions agree very well with the results of Pisa group jOj. 
Like in the n- 3 H case, MT I-III was found to be successful in describing differ- 
ential cross sections up to d+p+p threshold (see Fig. |2J). 

3.3 p+ 3 H scattering at very low energies 

He continuum is the most complex 4N system, since its spectrum contains nu- 
merous resonances (see Fig.^). Calculations of p+ 3 H scattering are furthermore 
complicated by the existence of the first J n = + excitation of 4 He in its thresh- 
olds vicinity. This resonance located at Er ~ 0.4 MeV above p+ 3 H threshold 
covers with its width r ~ 0.5 MeV almost the entire region below n+ 3 He. In 
order to separate n+ 3 He and p+ 3 H channels one should properly treat Coulomb 
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Figure 3. Calculated n+ 3 H total cross sections compared with experimental data 



interaction, the task is being furthermore burden since both thresholds are de- 
scribed by the same isospin quantum numbers. 

When ignoring Coulomb interaction, as was a case in the most of the calcu- 
lations performed until now, n+ 3 He and p+ 3 H thresholds coincide. + resonant 
state moves below the joint threshold and becomes a bound state for all NN 
potentials models we use. Former fact is reflected in low energy scattering ob- 
servables: on one hand N+NNN scattering length in + state is found positive, 
on the other hand excitation function decreases smoothly with incident particles 
energy and does not show any resonant behavior. 

By properly taking Coulomb interaction into account, thus separating n+ 3 He 
and p+ 3 H thresholds, we have placed the 4 He virtual state in between. However 
unlike in the other 4N systems, MT I-III predictions for p- 3 H singlet (J n = + ) 
scattering length as well as the excitation function - ^ (E) | g=1 - are in dis- 
agreement with experimental data. In this model 4 He virtual state is obtained 
too close to the p- 3 H threshold and therefore has very small width. Neverthe- 
less out of resonance region the excitation function approaches experimental 
data points, which is not a case in the calculations without Coulomb interaction 
considered. Our calculations with realistic potentials are still limited in PWB. 
Nevertheless, they provide very promising results displayed in Fig. 0] Pure 2NF 
models predict too large singlet scattering length, thus placing the virtual state 
too far from threshold. On the other hand, by implementing UIX-3NF in con- 
junction with Av.18 NN model one obtains singlet scattering length as well as 
the excitation function ^(-E , )| e _ 120o in agreement with experimental data |12j . 
A detailed analysis of these calculations is in progress ^Ql E] as well as their 
extension above the n-3He threshold. 
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